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I. INTRODUCTION 



Newly observed hadrons, in particular in the heavy flavor sector, have raised strong interest in treating and in- 
terpreting these states as hadronic bound states or, as commonly named, hadronic molecules (for overview see e.g. 
Ref. One of the main reasons to treat these observed states as molecules is that their masses are close to the 

thresholds of corresponding hadronic pairs. Canonical examples are the scalar D*g(2317) and axial Dsi(2460) mesons. 
As stressed in the scalar £'*q(2317) and axial Dsi{24:60) mesons could be candidates for a scalar DK and a axial 
D*K molecule because of a relatively small binding energy of ^ 50 MeV. 

In a series of papers 0-0 we developed the quantum field approach based on phcnomenological Lagrangians for 
the treatment of hadrons as bound states of lighter hadrons — hadronic molecules. In particular, we considered the 
strong, radiative and leptonic decays of £'*q(2317) and Dsi{24:6Q) mesons [3]-[5]. A new feature related to the DK 
and D*K molecular structure of the D*q and Dsi mesons was that the presence of u{d) quarks in the Z?^*-' and K 
mesons gives rise to direct strong isospin-violating transitions D*q — > Dgi:^ and Dgi —> D*Tr^ in addition to the decay 
mechanism induced by 77 — tt*' mixing, as considered previously. We showed that the direct transition compares with 
and even dominates over the 77 — tt" mixing transition in the isospin-violating decays of D*g and Dgi mesons. The 
composite (molecular) structure of the D*q and Dgi mesons was defined by the compositeness condition Z — 
(see also Refs. [ll| and [l]-[3l)- This condition implies that the renormalization constant of the hadron wave function 
is set equal to zero or that the hadron exists as a bound state of its constituents. The compositeness condition was 
originally applied to the study of the deuteron as a bound state of proton and neutron Q . Then it was extensively 
used in low-energy hadron phenomenology as the master equation for the treatment of mesons and baryons as bound 
states of light and heavy constituent quarks (see e.g. Refs. [1, [l^). We found that our theoretical framework gives 
numerical results which are consistent with both the experimental data and previous theoretical approaches. 

In this paper we extend our formalism 0-0 to the scalar i3*Q(5725) and axial Bsi{5778) bottom-strange mesons 
which have been considered before theoretically [l3|-[2^ as the bottom partners of the charm mesons D*o and Dgi- 
In Ref. masses, strong and radiative decays of B*q and Bgi states have been analyzed using heavy hadron chiral 
perturbation theory. In particular, the mass spectrum, strong and radiative decays widths have been calculated 
using quark models, QCD sum rules, effective field approaches based on chiral and heavy quark symmetry including 
coupled-channel unitarity. 

We assume that the i3*g(5725) and i3si(5778) are bound states of BK and B*K mesons, respectively. We adopt that 
the isospin, spin and parity quantum numbers of the B*q {B*q) and Bgi (Bgi) are I{J^) — 0(0+) and I{J^) — 0(1+), 
while for their masses we take the values mB^^ ~ 5725 MeV and ms^i = 5778 MeV (the central values predicted in 
Refs. 0, |2l|)- Note, that different approaches result in the following ranges for the B*q and Bgi masses: itlb*^ — 
5627 - 5841 MeV and ms^i = 5660 ~ 5859 MeV. Using a phcnomenological Lagrangian approach we analyze the 
strong B*^ -Bs7r°, Bgi — > B*tt'-^ and the radiative B*^ B*j, Bsi Bs-f, Bsi — > -8*7, Bsi — > -B*o7 decays. We 
give predictions for the decay properties: effective couplings and decay widths. 

In the present manuscript we proceed as follows. First, in Sec.|lTl we outline our framework. We discuss the effective 
mesonic Lagrangian for the treatment of the B*q and Bgi mesons as BK and B*K bound states, respectively. In 
Section III we consider the matrix elements describing the strong and radiative decays of the B*^ and Bsi mesons. We 
discuss our numerical results and perform a comparison with other theoretical approaches. In Section IV we present 
a short summary of our results. 



In this section we discuss the formalism for the study of the i3*g(5725) and Bsi (5778) mesons as hadronic molecules, 
represented by a BK and B*K bound state, respectively. We adopt that their quantum numbers (isospin, spin and 
parity) are I{J^) = 0(0+) for i?*Q (-8*9) and I{J^) = 0(1+) for Bgi {Bsi). For their masses we take the values 
tob*q = 5725 MeV and niB^i — 5778 MeV (the central values predicted in Refs. 0, Hlj). Our framework is based on 
an effective interaction Lagrangian describing the couplings of the i?*Q and Bsi mesons to their constituents; 



II. THEORETICAL FRAMEWORK 



A. Molecular structure of the Bso(5725) and Bsi(5778) mesons 




(lb) 



(la) 
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where the doublets of B, B* and K mesons are defined as 



B 



B+ 
B° 



B* 



B*+ 



K 



K- 



(2) 



The summation over isospin indices is understood and the symbol T refers to the transpose of the doublets B and B* 



The molecular structure of the B*^ and Bsi states is (we do not consider isospin mixing): | 
\Bsi) = \B*~^K^) + \B*^K°). In Eq. (P) we introduced the kinematical parameters w, 



B 



sO/ 



\B+K- 



mi/{mi + rrij), where 



i^B, TUB* and rriK are the masses of B, B* and K mesons. The correlation functions with M = B 



sO 



or Bsi characterize the finite size of the B*q and Bgi mesons as BK and B*K bound states and depend on the 
relative Jacobi coordinate y with, in addition, x being the center of mass (CM) coordinate. Note, that the local 
limit corresponds to the substitution of $m by the Dirac delta-function: $M(y^) S'^iu)- A basic requirement for 
the choice of an explicit form of the correlation function is that its Fourier transform vanishes sufficiently fast in 
the ultraviolet region of Euclidean space to render the Feynman diagrams ultraviolet finite. We adopt the Gaussian 
form, ^m{Pe/-^'ai) ^ 6xp(~pI;/Am) ' ^'^^ ^^'^ Fourier transform of vertex function, where is the Euclidean Jacobi 
momentum. Here Ab»^ is a size parameter, which parametrizes the distribution of B and K mesons inside the i?*Q 
molecule, while Ab^i is the size parameter for the B^i molecule. For simplicity we will use a universal scale parameter 
Am = ^B^g = Ab^j. 

The coupling constants gB*„ and gs^i are determined by the compositeness condition [1, Q, which implies that the 
renormalization constant of the hadron wave function is set equal to zero: 



Zb:, = l~S'B.^(m|.J = 0, 



(3a) 
(3b) 



Here, I]^.^(m^*^) = g^^ H^.^(m^*^) is the derivative of the B*q meson mass operator. In the case of the Bgi meson 
we have ^('^b i) ~ 3b ^'b i(™b i)' which is the derivative of the transverse part of its mass operator 
conventionally split into transverse 'Sb,! and longitudinal parts as: 



2 ^B, 



(4) 



where 



5± = • 



(5) 



The mass operators of the B*^ and Bsi mesons are described by the diagrams of Fig.[TJa) and[ljb), respectively. 
Following Eqs. (|3a|) and (j3b|) the coupling constants g^, and g^^^ can be expressed in the form: 



9l. 



(47rAM)2 
2 

(47rAM)2 



dx 





1 oo 



dx 



da a PQ{a, x) 
{l + af 

daa Pi(a,x) 
(l + a)3 



dzo 



2^^.(1 + a) dzi 



(6a) 
(6b) 



where 



PQ{a,x) 
Pi{a,x) 



a^x{l — x) + w^^^ax + i«^^a(l — x) , 
a^xil — x^+w^ ax + ail — x) . 



fi%ax + /i^a(l — x) 

li^,ax + /i^a(l — x) — 

rriM 
Am ' 



Poia,x) 2 
1 + a 
Pi{a,x) 
1 ^ 



(7) 



The above expressions are valid for any functional form of the correlation function $ 



M- 
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Note that the compositeness condition of the type (|3ap . (j3bp was originally applied to the study of the deuteron as a 
bound state of proton and neutron [Sj . Then this condition was extensively used in low-energy hadron phenomenology 
as the master equation for the treatment of mesons and baryons as bound states of light and heavy constituent 
quarks ^,^0*1. In Refs. [IH and 3-0 this condition was used in the application to hadronic molecule configurations 
of light and heavy mesons. 

B. Effective Lagrangian for strong and radiative decays of B*q and Bsi mesons 

In Ref. f3l,l5|, in the analysis of the strong isospin-violating decays D*q I?s7r° and Dgi -D*7r° in the molecular 
approach, we showed the existence of two possible dynamical mechanisms. This included the so-called "direct" 
mechanism with 7r°-meson emission from the D D* and K K* transitions and the "indirect" mechanism where 
a TT*^ meson is produced via rj — Tv^ mixing. The mixing is due to the mass term of pseudoscalar mesons in the leading- 
order 0{p^) Lagrangian of chiral perturbation theory (ChPT) [2^, [sO]. Note, that the second mechanism based on 
ry— TT*^ mixing was mainly considered before in the literature. Originally it was initiated by the analysis based on the use 
of chiral Lagrangians [sO, [Sll, [s^l where the leading-order, tree-level D^qDsTt^ (DsiDItt^) coupHng can be generated 
only by virtual 77-meson emission. In Ref. Q we showed that the "direct" and "mixing" mechanisms can be combined 
together in the form of an effective coupling of tt'^ to the mesonic pairs DD* or KK* with modified flavor structure. 
This modification occurs after the diagonalization of the mesonic mass term involving tt'^ and 77 meson fields ^29j (see 
details in [^). In particular, instead of the r^n^ coupling to DD* or KK* we have 7r° (ra cose + kI sine), where 
K = l/-\/3 or -\/3 is the corresponding flavor-algebra factor for the DD* or KK* coupling, respectively. The ~ rj 
mixing angle e is fixed as [1^: 

tan2e=— ^ ~ 0.02 , m = -(to„ + m^) , (8) 

2 iris ~ m 2 

where to„, nid, rrig are the current quark masses. 

Below, in Eq. fTTjl . we display the explicit form of the corresponding interaction Lagrangian. The lowest-order 
diagrams which contribute to the matrix elements of the strong isospin-violating decays B*q Bgir^ and Bsi S*7r'^ 
are shown in Figs. 2 and 3. Note, that in the isospin limit (m„ = m^), the 77 — 7r° mixing angle vanishes and the 
masses of the virtual i?*^*-* and if^*^ mesons in the loops are degenerate. As a result the pairs of diagrams related to 
Figs. 2(a) and 2(b), Figs. 2(c) and 2(d), Figs. 3(a) and 3(b), Figs. 3(c) and 3(d) compensate each other. Therefore, in 
the calculation of the diagrams of Figs. 2 and 3 we go beyond the isospin limit and use the physical meson masses. 

The diagrams contributing to the radiative decays B*q B*j, Bgi Bsj, Bgi B*^ and B^i B*qJ are 
shown in Figs. 4, 5, 6 and 7. The diagrams of Figs. 4(a), 4(b), 5(a) and 5(b) are generated by the direct coupling 
of the charged K^ and B^{B*^) mesons to the electromagnetic field after gauging the free Lagrangians related to 
these mesons. The diagrams of Figs.6(a)-6(d), 7(a) and 7(b) are generated by the coupling of a corresponding pair 
of vector and pseudoscalar mesons to the photon. The diagrams of Figs. 4(c), 5(c), and 6(e) are generated after 
gauging the nonlocal strong Lagrangians ^ describing the coupling of the B*q and Bsi mesons to their constituents. 
The diagrams of Figs. 4(d) and 5(d) arise after gauging the strong B*B~^K~ and BsB*^K~ interaction Lagrangian 
containing derivatives acting on the pseudoscalar fields. Details of how to generate the effective couplings of the 
involved mesons to the electromagnetic field will be discussed later. 

After the preliminary discussion of the relevant diagrams, we are now in the position to write down the full effective 
Lagrangian for the study of the strong and radiative decays of the B*g and Bgi mesons formulated in terms of mesonic 
degrees of freedom and of photons. We follow the procedure discussed in detail in Refs. 3-0], where we considered 
the I?*Q and Dgi meson decay properties. First, we write the Lagrangian C, which includes the free mesonic parts 
£freo and the strong interaction parts £str: 

C{x) ^ £froc(a;) + >Cstr(a;) , (9) 

with 

£free(x) = -B: + ix)iD + ml,jB:,ix) + P,i.^(a:)(g'-[n + mlj - d^^d^BsiA^) 

- l^[x){U+ml)7:{x) + ^-^[n\x)f- P\x){U + ml)P{x) + ^ 5pP\x)P\x) 

P=K,B,B^ P=K,B 

+ E V^{x){g^^''p+ml]~d'^dnK{x)~ SvV^{x)V'^^{x), (10) 

V=K*,B*,B* V=K',B' 
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Cstrix) - iMl^B;^{x)7^Bix)^^ Bix) + K;\x) 7:k{x) td K{x) 



+ g,^,,,B;\x)K{x)iTB^,{x)+g,^,,,K;\x)B{x)iTmx)+9s>,.B:l{x)B^{x)iTK{x) 



BI^^-{x)BI\x)K:{x)+BII{x)K*^{x)B:^^{x)+K*^'^{x)B:1{x)BI\x) 
e'^'^"^B:^,{x)B:l{x)K{x)+CB:,{x)+CBAx) + H.c. , (11) 



where summation over isospin indices is understood, □ — d'^d^ and A d B = AdB — BdA. Here, 7? = (tti, 7r2, tts) is 
the triplet of pions, ttb = ttiTi +772X2 + 773 (ts cose + / sine/ Vs), TTi^ = ttiti +7r2T2 + 7r3(r3 cose + /sine \/3), /?'*' and 
are the doublets of pseudoscalar (vector) mesons, Bf and B* ± are the pseudoscalar and vector bottom-strange 
mesons, respectively, T/*^" — d'^'^V*'^ — d^V* ^ is the stress tensor of the vector meson field. 

In our convention the isospin-symmetric meson masses of the isomultiplets m^^mp^my are identified with the 
masses of the charged partners. The quantities 5m are the isospin-breaking parameters which are fixed by the 
difference of masses squared of the charged and neutral members of the isomultiplets as: 5m = ™m± — ™a/o ^^^^ 
niMO = TOjgo . The set of mesonic masses is taken from data [s^ . From Eq. ((TT|) it is evident that the couplings of tt*^ 
to the B* B and K*K mesonic pairs contain two terms — the "dominant" coupling (proportional to cose) and the 
"suppressed" coupling (proportional to sine). This means that the first coupling survives in the isospin limit, while 
the second one vanishes. 

The free meson propagators are given by the standard expressions 

^DM{x-y) = {Q\TM{x)M^{ym j ^^^e^^^-v) D,,{k) (12) 
for the scalar (pseudoscalar) fields, where DM{k) — ~ k"^ ~ ie)^^ and 

(2^i 

for the vector (axial) fields, where D'^^,{k) — (— g'"' + fc^fc'^/m^j,) (to|,/* ~ — ie)^^ . The choice for the strong meson 
couplings of the Lagrangian (jlip will be discussed in Sec. III. 

The electromagnetic field is included in the Lagrangian ^ using minimal substitution i.e. each derivative acting 
on a charged meson field is replaced by the covariant one: a/^M(*)± ^ (9^ =Fieyl^)M(*)±. Note, that the strong 
interaction Lagrangians Cb*„ and £_b,i should also be modified in order to restore electromagnetic gauge invariance. 
It proceeds in a way as suggested in Rcf. [34] and is extensively used in Refs. [l3]- In particular, each charged 
constituent meson field //^ (i.e. /?(*='=) and K^'^'>) in and Cb,i is multiplied by the gauge field exponential (see 
further details in [i,[i,[i3]): 

H^iy) e^"''^y''''^^H^{y) (14) 

where 

X 

I{.x,y,P) = J dz^A^(z). (15) 



z/?g4^-2/)-(0|TA/*^(.T)A/*^t(y)|o) = / ^^e-'^^--^) ^^^.(fc) (13) 



For the derivative of I{x,y,P) we use the path-independent prescription suggested in 3J] which in turn states 
that the derivative of I{x,y,P) does not depend on the path P originally used in the definition. The nonmini- 
mal substitution (|14p is therefore completely equivalent to the minimal prescription. Expanding the exponential 
gT«e^(a,a:,P)^±^y') powcrs of the electromagnetic field and keeping linear terms like the four-particle coupling 
B*qB^K~j and BsiB*^K^j we generate the diagrams of Figs. 4(c) and 5(c), 6(c), 6(d), respectively. 

Finally, we specify the electromagnetic Lagrangian describing the coupling of vector and pseudoscalar mesons to 
the photon: 

CvP^{x) = ^F^,ix)e'^'''^^(^g^,^^^^K: + {x)K-{x)+g^^,,^„^K:^^ix)K°ix) 

+ g,,^,^_^B*Jix)B-ix)+g^,,^,^B*^°p{x)B°ix)) + H.c. (16) 
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Here the couplings Qyp^ can be extracted from the corresponding decay widths V — s- P + 7. Presently we only have 
information about the decay width of K* mesons. Using the expressions for the K* K + j decay widths 



and the data (central values) for r{K*^ K^j) = 50.29 keV and T{K*° -> K°'y) = 116.19 keV we deduce the 
coupling constants ^ = 0.836 GeV^^ and g^.o = —1.267 GeV^^. Note, that in the nonrelativistic SU(3) 
quark model the couplings ^± and 3^,0 ko proportional to the sum of the charges of the constituent 

quarks: .g^,,± ^± ''^ (e„ + Cg) = 1/3 and ff^.o^^^o.^ ipd + e^) = —2/3. This is the reason why the coupling of 
the neutral kaons is defined (by convention) with a negative sign. Also, the prediction of the nonrelativistic quark 
model for the ratio (7^,0 ^„ /3^,,± ^± = —2 is violated by relativistic corrections. For D mesons the corresponding 
ratio is in precise agreement with data. In particular, taking the experimental values of g^^^ ~ 0.5 GeV~^ and 
5£,.±__D±.^ — 2-0 GeV~^ (see the discussion in Ref. [Ill) we get 

^^•"-""^ = + = 4 (18) 

9o'±^o±-, + 

Therefore, one can expect that for bottom mesons the naive quark model should prediction should also be sufficient: 

^B'O^BOy Sh + ed 



= -2. (19) 



For our numerical estimates we will use typical values with g± — 0.5 GeV ^ and g^^^ = — 1 GeV ^. 

These couplings correspond to the full width of B*^ equal to 0.23 keV and and of equal to 0.91 keV (as usual 
we suppose that B* B^ is the dominant mode for the B* mesons). 

III. STRONG AND RADIATIVE DECAYS OF THE B*o AND Bsi MESONS 

A. Matrix elements and decay widths 

The matrix elements describing the strong B^q Bs7!'''\ Bsi B*^'^ and radiative B*q — > -8*7, B^i Bgj decays 
are defined as follows 

M{B:,{p)^Bs{p')n''{q)) = Gbi,b^.. (20a) 

M{B,,{p) ^ B:{p')n"{q)) = e^ip)e:{p') {g^'' Gb.,B'. ~ v'^v'' Fb^.b-.) . (20b) 



and 



M{B:^{p) ^ i?:(p')7(g)) = e el{q)el{p') (g^^p'q - p'^q,) Gb;„b;, , (21a) 
M{Bs,{p) ^ B,{p')^{q)) - e e^{p)e:{q) (g^^ pq - q^p^ Gb^.b^j , (21b) 

MiB,,{p) B^ip'hiq)) = ee"™'"^e"(p)e*^(p') £;(<?) <Z. (gb^.b*^ g^„5amM 



+ Fb^iB^-Y gp,nPmqa + Hs.iB^-f gam Puq^i j , (21c) 

M{B,,ip) ^ B:„{p'h{q)) - e e^{p)et{q) p^ qp Gs.1B.07 > (21d) 

where v = p/mB^^ and v' — p' /iriB* are the four-velocities of the Bsi and B* mesons, Gb'^b^tt, G{F)B^iB^n, ^B'^B^-y 
and Gb^iBsjj G{F, H)B,iB^y and Gb^iB^o^j are the corresponding effective coupling constants. The coherent sum of 
all the diagrams in Figs. 4-7, contributing to the radiative decays of B*q and Bgi mesons, is gauge invariant, while 
the contribution of each diagram is definitely not gauge invariant. As done in Ref. for convenience we split each 
individual diagram into a gauge-invariant piece and a remainder, which is noninvariant. One can prove that the sum 
of the noninvariant terms vanishes due to gauge invariance. In the following discussion of the numerical results we 
will only deal with the gauge-invariant contribution of the separate diagrams of Figs. 4-7. In Appendix A we present 
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the calculational technique for determining the effective couplings entering in the matrix elements of the strong and 
radiative transitions of B*q and i?si mesons. 

Using Eqs. ^ and ^ the strong B*o B^n", B,i -> B*t:° and radiative B*q -> B*7, B^i 5^7, B^i B*-/, 
Bsi — > B*Q^ decay widths are calculated according to the expressions: 



V{B:,^B,^^) = ^^rf^^.i, (22a) 



sO 



and 



nBsi ^ = + ^ (gs.,b> - - 1) j | , (22b) 



T{B:^ ^ S:7) = a G|.^s.^ , (23a) 
T{Bsi^B,-i) = ^ G|^,s.^ P^^a , (23b) 



T{Bsi ^ S:7) = I P^'a { (gs.iB.-, + ^s^iSj^) + ^ (gs^.b.^ + } . (23c) 

r(i?,i^B:„7) = ^Gl^^s.^^P^^, (23d) 

where w = vv' = (^s^i + "^b* ~ m^o)/ (2mB^i'rnB* ) and Ptti, Pyi are the corresponding three- momenta of the decay 
products. 

Note, the contribution of the effective coupling constant Fb^iB-it to the Bgi B*7t'^ decay width is strongly 
suppressed. This is because the contribution of the matrix element with Fb^^bitt is proportional to the suppressed 
factor w^ — 1^4x10^"^ with w ~ 1. Therefore, we have 



and 



r(P;o -> P.ttO) P^i\mB^J \Gb'b. 



(25) 



B. Numerical results 

First, we discuss the choice for the strong coupling constants in the Lagrangian £str In Refs. J^J^ we used the 
set of strong coupling constants g^^o^L {9d'd^^ 9d,d'k = 9d,k-d and g^.^^ = go'o'K* ) defined in the charm sector, 
where index L denotes a light meson, while Di and D2 are the respective charm states. The coupling g^,^^ — 17.9 
was deduced using data for the corresponding strong decay width [36[. The coupling constants g^^ ^.^ and ff^.^^ have 
been estimated using two different variants of the QCD sum rule approach discussed in Refs. [sjlll], where similar 
results have been obtained. Both Refs. [sl, HI] point to a strong suppression of these constants in comparison to the 
coupling 5^3. o„ ■ An updated analysis for g^,^^ in the context of a QCD sum rule approach gives a result close to data 
- g^^o^ = 14± 1.5 [39[. We used the predictions of Ref. |37|]: 5^.^, ^ = 2.02 and g^,^j^ = 1.84. For the unknown 
parameters g^^^^,^ and g^,^,^, we used the approximative relations 5^^^.^ 9o,d'k and g^.^'K' - 9otDK^ vfhich 
can be explained phcnomcnologically: the first relation - by the universality of the coupling of the Ds meson to D* K 
and K*D mesonic pairs (it is based on exact SU(4) flavor symmetry and we do not expect a substantial violation of 
this relation due to breaking of the SU(4) symmetry) and the second relation - by the universality of the coupling of 
the D* meson to two pseudoscalars and two vectors (like for ptttt and ppp couplings: ^pTrir — gppp — 6, and also for 
J/^DD and J^D*D*: gj/^oD — gj/^D'D' — 8 O)- For consistency, in the present manuscript we use the set of 
charmed hadronic couplings predicted by QCD sum rules (central values) with : 

9o'n. = 14 , go.n-.- = 9n^,^.n = 2.02 , g^,^^ - g^,^,^,, = 1.84 . (26) 

In particular, instead of g^,,^^ — 17.9 we use ^^.^ir ~ Such a modification does not change the numerical results 
of Refs. [3, Q, because the contribution of the corresponding diagrams containing the coupling g^,^^ is strongly 
suppressed. A similar picture we also have in the bottom sector (see discussion below). 
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In order to evaluate the corresponding bottom couplings gB_^g^^ we use the arguments of heavy hadron chiral 
perturbation theory (HHChPT) [4l|, |43| which relates the bottom and charmed couplings containing the same light 
meson: 

9biB2L 9diD2L l^'J 

where ms and m^i are the masses of the bottom and charm mesons identified e.g. with the masses of mg± = 5.279 GeV 
and TO£i± = 1.8693 GeV. In other words, to get the set of bottom meson couplings used in the strong Lagrangian (fTTj) 
we rescale the corresponding charm meson couplings by a factor mg± /mjj± ~ 2.82. In particular, we have: 

9b-b. = 39.5 , - 9b,k'b = 5.70 , g^.^^ ~ 9b;b*k' = 5.19 . (28) 

For the coupling gs'B'K we have no further input. From a dimensional analysis it should be of order 1 GeV^^. The 
diagrams 6(e) and 6(f), however, where this coupling enters, are strongly suppressed, and we therefore do not need 
to have precise knowledge of this constant. 

The coupling g^,^^ = 4.61 is fixed from data on the K*Kn decay width [s^. Finally, we fix the couplings g^, 
and gg^_^ , which are given by Eqs. ([6a|) and (|6b|) in terms of the adjustable vertex function. Using the Gaussian vertex 
function we obtain the result that these couplings are quite stable with respect to a variation of the scale parameter 
Am- In particular, varying Am from 1 to 2 GeV, we get a range of values for g„, from 27.17 to 23.21 GeV and for 
g^ ^ from 25.64 to 22.14 GeV. Note that our predictions for these couplings are in agreement with the results of other 
theoretical approaches: g^, = g^^^ = 19.6±5.7 GeV [Hght-cone QCD sum rules approach ^] and g^, = 23.572 GeV 

and g^^^ = 23.442 GeV [effective chiral approach 0,123]. 

Now we present the numerical results. First, we discuss the results for strong decays. Here the main contribution 
to the decay width comes, as expected, from the diagrams of Figs. 2(a), (b) and 3(a), (b). On the other hand, the 
contribution of the direct mechanism is comparable to the one of the indirect mechanism (i.e. due to ry — tt'^ mixing). 
The contribution of the diagrams in Figs. 2(c) and (d) is of order 0.1% of the total contribution to the Gs^gBsTr 
couphng and the contribution of Figs. 3(c) and (d) is of order 1.6% of the total contribution to the Gb^iB^u coupling. 
Therefore, the couplings Gb'^BsTt and Gb^iB*-w are not sensitive to a variation of the couplings gg.^^, 9b* bk 
g^ j3*/f ■ They are only sensitive to the values of the couplings g^ and g^^g,^, ■ The leading-order contributions 
from the diagrams in Figs. 2(a), (b) and 3(a), (b) to the quantities Gb'^b^tt and Gb^iB^tt in terms of the couplings 
Sb^k-b and g^.^.^^. are given by 

Gb:,b,. = 65.9 gg^^,g MeV , Gs,,b-x = 72.3 gg,g,^, MeV (29) 

for a typical value of the dimensional parameter Am — 1 GeV. Then, using the specific values of gg ^,g = 5.70 and 
9b'b*k' = 5.19 we get the following predictions for the effective couplings and decay widths: 

Gb'^b.tt = 375.7(375.6) MeV , Gs^iS^ = 381.1(376.3) MeV (30) 

and 

r(S;o ^ B,7r°) = 55.2(55.1) keV, r(S,i ^ B>°) = 57.0(55.6) keV . (31) 

In the brackets we indicate the results of the leading diagrams of Figs. 2(a), (b) and 3(a), (b). 

The strong decay couplings Gb'^b^-w, Gb^iB'-k and the decay widths r(i3*Q — > -B^tt), T{Bsi — > -B*7r) are practically 
degenerate, which can also be explained by heavy quark symmetry (HQS), which is a good symmetry for the heavy- 
light mesons with a bottom quark. Because of the infinitely heavy mass of the bottom quark the spins of the h 
antiquark and the s quark decouple (spin symmetry), and, therefore, the properties of i3*o and B^i mesons become 
similar. This issue was also discussed in Refs. [3l|, [S^l in the context of the charm partners (I?*o and Dsi). In 
particular, it was shown that the corresponding coupling constants and widths are degenerate in the heavy quark 
limit. In our previous papers we also reproduced the same result. Moreover, for finite charmed meson masses the 
decay characteristics are nearly degenerate: 

Gdi,,d,^ = 145.4 MeV , Gd^^d^^ = 160.2 MeV (32) 

and 

r(D;o -> DsTr°) = 46.7 keV , r(D^i D>°) = 50.1 keV . (33) 
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The next point is a check of the flavor content of HQS. In our molecular approach the leading contributions to the 
couplings Gn'gHsTr f^nd GHsiH;Tn where H — D or B, are defined by the pairs of diagrams in Figs.2(a,b) and 3(a,b), 
respectively: 

Gh'^^h^tt = gH^^9K*K„ 9„^HK- ^KK'H , (34a) 
Gh,iH;-k = gH.^aK'K^aH'H'K* ^KK'H' , (34b) 

where Ikk'H and Ikk'H* are the structure integrals which are of order ©(l/mg) in the inverse heavy quark mass 
expansion - 1/toq. The hadronic couplings scale as: 

VJo ^ ' 9H,^ ^ O(toq) , g^,^^ - 0(1) , g„^„^, 0{mQ) , g„,„,^, ^ 0{mQ) . (35) 

Therefore, the strong couplings Gh^^Hstt and Gh^iH'tt scale as: Gh^^jHstt ^ 0{raQ) and Gh^iH^tt ^ 0(mQ). For this 
scaling behavior the following relations result: 

Gb'B.w nib GbsiB^tt nib 

jr^ , 7^ ^ ■ (36) 

With our results [see Eqs. (|30p and ((32)l ] we conclude that the constraints ((36| are fulfilled very well. 

In Table [J we present our results for the decay widths r(_B*Q BsTt^) and r(i?si ^ BItt^) including a variation 
of the scale parameter Km from 1 to 2 GeV (an increase of Km leads to an increase of the widths) and compare 
them with known theoretica,! predictions [20l . [2ll Wf\ . Our results for the decay widths are larger in comparison 
to previous approaches [H, [2l|, [l^l due to inclusion of the direct isospin-violating transitions i3*Q Bgir^ and 
Bsi ^ 

Now we turn to the discussion of the radiative decays B*q — s- B*j and Bsi — > Bg^. The main contribution to the 
decay characteristics of the B*q B*'-f and Bgi Bgj decays comes, as expected, from the diagrams of Figs. 4(a) 
and 5(a). Our results for the effective coupling constants and decay widths for a typical value of Am = 1 GeV are: 

Gb*„s-7 = 0.122 GeV-S Gb.iB,^ = 0.115 GeV"^ (37) 

and 

T{B;o^B;j) = 3.07 keV, T{Bsi B,"f) = 2.01 keV. (38) 

In Table [TT] we summarize our results for the radiative decay widths including a variation of the scale parameter Am 
from 1 to 2 GeV (an increase of Kbm leads to a larger value for the width). In comparison, we also display the 
predictions of other theoretical approaches [H, [H, . The lower limit of the QCD sum rule results [1^ is consistent 
with our predictions. In our opinion the predictions of [H, are overestimated. In particular, applying HQS (spin 
symmetry) we can relate the corresponding radiative coupling constants of the same flavor as 

G_D»|.j_Dj7 = Gd^iD^j , ^B^gB'-i — Gb^iBs'j ■ (39) 

The same relation was derived previously in HHChPT 32] in the charm sector. Now we perform the same exercise as 
done for the strong couplings in order to relate the radiative couplings of different flavors. The leading contributions 
to the couplings Gn'^H^-y and GH^iH^-y, are defined by the pairs of diagrams in Figs. 4(a) and 5(a), respectively: 

Gfl'gH^i = gn'g 9kk-, gH*HK Ikkh , (40a) 
GH,iHs-t = 9h,i 9kk-, gH.H'K Ikkh- , (40b) 

where Ikkh and Ikkh* are the structure integrals which are of order 0{l/mQ). The couplings g^K-,-! 9h*hk ^^'^ 
9hsH-k scale as: g^^^ ~ 0(1) , ~ O(mQ) and g„^„,j, ~ 0{mQ). Therefore, the radiative couplings Gui^Ht-t 

and Gh^iHsj are insensitive to the flavor of the heavy quark/meson: Gh'^^h^j ~ G{1) and Gn^iH-y ^ 0{l). Finally, 
the radiative couplings obey both spin and flavor symmetry in the heavy quark limit and we arrive at the constraint: 

Go-gD'y — Gd^iD^j = Gb'^B^j — Gb^iB^-j ■ (41) 

Recalling the results for the radiative decay constants of charmed mesons [1] 

Gdi,,di-i = 0.093 GeV^S GD,^D,■y = 0.106 GeV~^ (42) 
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we conclude that relation (|4ip is approximately fulfilled by our results, which are obtained for finite physical masses of 
the heavy mesons. The constraint ([1T|) can also be used to deduce relations between the corresponding decay widths: 



and 



TjDsi ^ Dsl) 



TjBsi ^ B,7) 





3.80 , (43a) 

(43b) 

3.74 , (43c) 

(43d) 



Relation (|43ap is confirmed by the full analysis done in our framework and in other theoretical approaches (see 
compilation of the results in Refs. 3, 5]). E.g. the relation explains why most of the approaches predict that the 
decay width T{Dsi Dg^) is approximately 3 — 5 times larger than V{D*q D*^). The other relations (j43bp - (|43d[) 
help to give predictions for the decay widths of the bottom partners. In particular, we arrive at the conclusion that 
our full predictions for r(i3*Q — > i3*7) and T{Bsi — > Bs^) ~ a few keV - are well justified. 

Now we discuss the ratios of radiative and strong decay modes. For both systems, i3*o and Bji, we predict small 
ratios: 

0.05, 

0.03 . (44a) 







^Btj 




BsTT 






-^Bsj 


Bsi - 


-.B^TT 



Note that similar results we also obtained in the charm sector: Rd*^ = 0.01 and Rd^i = 0.05. The predicted ratio 
Ro^g is consistent with the present experimental limit Rd^^ < 0.059, while the ratio Rdsi is smaller than the result 
quoted by the Particle Data Group Rd^i — 0.38 ±0.05 [33j. We consider this experimental result as preliminary, since 
it is not clear why the ratio for the axial state D^i is much larger than for the scalar state D*q. At this point more 
precise data on the strong and radiative decays of the D^i meson and their ratio would be very helpful. 

Finally, we give the predictions for the other two radiative decay modes of Bgi mesons - Bsi — > B*^ and Bsi B*^^. 
These decay amplitudes are generated by the anomalous couplings of two vectors and one pseudoscalar and one can 
expect that they are suppressed. Moreover, there is an additional mechanism for their suppression. The leading 
diagrams in the process Bsi — > -8*7 are the ones of Figs. 6 (a) and (b). The separate contribution of the diagrams 
in Figs.6(c)-(f) to the corresponding decay width is of order 1%. The diagrams of Figs. 6(a) and Figs. 6(b) are 
subtracted from each other because of the opposite sign of the anomalous couplings 5^.±^j^± = 0.836 GeV~^ 
and g^.o ^0 = —1.267 GeV^^. In the case of the Bgi — > -S|g7 decay we only have two diagrams contributing 
to the matrix element. Again, due to the opposite sign of the anomalous couplings = 0.5 GeV~^ and 

9b*o bo ~ ~^ GeV^^ their total contribution is given by the difference of the individual contributions. Finally, as 
a full result we get values of 0.04 to 0.18 keV for the decay width of Bsi B*j including a variation of the model 
parameter A from 1 to 2 GeV. The result for the decay width of Bsi B*q'j is stable with respect to a variation of 
the model parameter A in the same range and is equal to 0.022 keV. We also present the result for r(Bsi -S*o7) in 
terms of the couplings , , and and „ „ : 

r{Bsi ^ b:,j) - [(5,.±^,±, +5«.o^,o,) X 1 GeV]2 0.088 keV. (45) 

Recently, these radiative decay widths have also been estimated using light-cone QCD sum rules [28j: r{Bsi 
B*j) = 0.3 - 6.1 keV and T{Bsi -> B^^-f) = 0.002 - 0.008 keV. The decay width T{Bsi B;^j) is definitely strongly 
supressed in both approaches. The decay width r{Bsi B*^) predicted in [28] is also smaller in comparison to the 
other two modes B*q i?*7 and Bsi — > -Bs7 (see Table HI)) . 
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IV. SUMMARY 

In this paper we studied the new bottom-strange mesons B*q{5725) and i?si(5778) in the hadronic molecule inter- 
pretation, i.e., we considered them as bound states of BK and B*K mesons using a phenomenological Lagrangian 
approach. Our approach is based on the compositeness condition with Z = (the renormalization constant of the 
molecular state equals zero). This condition is crucial for weakly bound states. The compositeness condition is a 
self-consistent tool to evaluate for a hadronic molecule, once the mass of the bound state is fixed, its coupling to 
the intermediate hadronic state, which in turn feeds the accessible final states. Once the further hadronic couplings 
to generate the final states are known, this approach is fairly model-independent in the results for the observable 
decay modes. Furthermore, the Lagrangian method implied by the compositeness condition allows a fully Lorentz 
and gauge-invariant treatment of the problem. 

We calculated the strong B*g BsTt'^, Bsi B*Tr° and radiative B*^ ^ B*j, Bsi Bsj, Bsi ^ B*j, 
Bsi SJq7 decays. A new impact of the BK and B*K molecular structures of the i?*o(5725) and i3si(5778) mesons 
is that the presence of u{d) quarks in the B^*^ and K meson loops gives rise to direct strong isospin-violating transitions 
B*g BsTt'^ and B^i — > B*tt'^ in addition to the decay mechanism induced by 77 — tt" mixing. We showed that the 
direct transition is comparable with the 77 — tt" mixing transition. As a consequence, the presence of the direct mode 
makes our predictions larger than the ones of previous approaches. In the case of the radiative decays B*^ — > B*^ 
and Bsi Bsj, our results are considerably smaller than in previous calculations. We also gave predictions for the 
anomalous decays Bsi B*j and Bsi — >■ Bsqj: their decay widths are suppressed in comparison to the leading 
radiative modes i3*Q — > Btj and Bsi Bs'j. 

In series of papers already studied in detail strong, electromagnetic and weak decays of D*q and Dsi 

states, as a consequence of their possible molecular structure. At present we do not think that the structure issue 
concerning the I?*q and Dsi mesons is settled yet, in particular because the mass and the narrowness of these states 
cannot be easily explained in the context of the standard cs picture [H, 0] ■ A further direct consequence and analogy 
of possible hadronic molecules in the Dg sector is also found in the Bg system. If the _B*g and Bsi mesons can possibly 
be experimentally established near the predicted mass values, then, because of their closeness to the BK and BK* 
thresholds, they are clear candidates for hadronic molecules. In the present approach we give clear predictions for 
the possible strong and radiative decay modes of -B*g and Bsi mesons accessible by experiment. Strong experimental 
deviations from our results would discard the molecular interpretation of these states. We hope that our results will 
be useful for future experiments, where _B*q and Bsi could possibly be detected. 
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APPENDIX A: EFFECTIVE COUPLINGS FOR STRONG AND RADIATIVE TRANSITIONS OF B*o 

AND Bsi MESONS. 

Here we discuss the calculational technique of the matrix elements of strong and radiative transitions of B*q and 
Bsi mesons. 

As an example for the calculation of diagrams presented in Figs. 1-7 we consider a generic loop integral containing 
a product of n virtual momenta fc^^ • • • fc^^, three meson propagators with masses mi, m2 and ma and the correlation 
function of the molecular state {B*q or Bsi meson) 

I (n n'^ = [ —S,(-p] ^Ati • • • ,-A1^ 

Vi-P„lP,Pj J ttH^ ' {mj - {k + p)^){ml - {k + p')^)iml ~ k^) ■ ^ ' 

The three main ingredients are 

• use of the Laplace transform of the vertex function 



which is useful to proceed with vertex functions of any functional form, 

• the a-transform of the denominator 

c 





• the differential representation of the numerator 

where i? is a linear combination of the external momenta. 

The calculation of the transition form factors amounts to a one-loop integration. Integration over the loop momentum 
is done analytically. One ends up with the integrals over a Feynman parameters which are not difficult to evaluate 
numerically. All calculations are done by using computer programs written in FORM [4^ and in FORTRAN for 
numerical evaluations. Also, for transparency we give explicit expressions for the leading contributions to the effective 
couplings defining the structure of the matrix elements of strong and radiative decays of B*q and Bsi mesons. 

1. Decay B*o ^ BsTt°. 

The leading contribution to the coupling constant Gb'^b^-k coming from the diagrams in Figs. 2(a) and 2(b) is 
defined by 

GB*gB,7T = 9B*^ gK'KTT gs^BK' ^B^gB.TT ■ (A2) 



Here Ib*qBs'k is the structure integral 

5 



i=l 

DC DO DC' 



Ibi,b,. = J J J Af exp(^i)iVi^ 





OC DO 'OO 





CO oo oo 



oo oo CJO 





OO oo 

//^ exp(A5)A^5, 



' 



where we introduce the set of notations: 

Ai = 1 + ai23 , A2 = A3 = A4 = 1 + ai2 , A5 = 1 + ai , 

Ai = -ain% - a2/x^. - as/^l- - I^b'WbkWkb + -r- i l^'^'aiKBOt^BK + IJ'\aiKBa2 + IJ?T,a3BKa2 

SU \ SU 

A2 = -ailJ,% - Q!2m|-- - IJ'%'^^WbkWkB + (ylJ'^'^^Q-XKBWBK + Ms,"lA'SQ:2 + ^Ji'iwBK0^2^ , 
A3 = -aill\, - a2HK - IJ^^I^WbkWkB + (^%*^a2BKWKB + IJ^Bs'^KBOil + l^la2BKai^ , 
Aa = a2lJi% - ^1%*WbkWkB + -T-IJ^%\(^\KB0i2BK , 

A = -aiii\, - ij?b*^wbkWkb + (^IJ'%'^^wbkWkb + h^^wkboii + nlwBKai^ , 

2 2 

iVi = iJ,% + i^K- 2/i|. + mI, +(4- ^J■%' + 2~ (mL - ' 

Ar2 = 1 + 
= 1, 



WM1M2 = ; I "123 = ai + a2 + as , ai2 = ai + 02 , aiMiMz = "i + wm^m^ , I^m = — ^ 

rriMi + A 
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2. Decay B^i ->■ S>°. 

The leading contributions to the coupling constants Gb*^b*tt and Fb*^b'tt coming from the diagrams in Figs. 3(a) 
and 3(b) are defined by 

GbIiBJtt = 9B,igK'K7rgB'B'K'lB,iB*Tv (A4) 

and 

Fbi^B'-k = 9Bsi9k'K7t9b*B'-K''IEsiB*-k ■ i-^^) 
Here Ib^iB*tt ^^'^ ^BsiB'-r structure integrals 

5 

tG _ V- TG,i 

^BsiB'n — / ,^B.,B'n ' 

i=l 

oo oo oo 

= ///^^^i^exp(Af)7Vf, 



oo oo oo 



oo oo oo 

iZb:. = j j j'-^^MA'i)N-, (A6) 



oo oo oo 





oo oo 



jG,4 
BsiB*7r 







and 



tF _ tF.i 

oo oo oc 

Ib\,bi^ = J J J Af ^""^^^ ^ ' ^ ^ 





cx:; oo DC' 

tF,^ _ III daida2 i aF\ -\tF 







where we introduce the set of notations: 



Ai = = -ai/u|f. - a2lJt\ - "sMs* - 1^'b,.,wb*kWkb' (^fjii^^otiKB'CtaB'K + A^0iiKB*a.2 + IJi'^a^B' Ka2^ 

A2 = A2 = -ail^B - Ct2l^K' - I^Bsi'^B'KWkB* + (^IJ'%^^aiKB*WB'K + IJL% ^Oi\K B* 0-2 + ^J^^^,WB^ KOi^ , 

A3 = -aill\, - Q.\li\ - IX^b^^Wb'-kWkB" + (^\^^a2B*KWKB' + li\^WKB*(y-l + /"^Q!2S*if Ctl^ , 

A^ = -Q!iAt| - ai/^K - I^Bsi'"^b*kWkb' + -^IJi%^^a\KB'a2B*K , 

A^ = -ai/xlf. - ij,%^^wb*kWkb* + \ H%^^wb'kWkb* + m|,wk-b*Q:i + /U^WB*ifQ!i ) , 
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ArG o 2 2 2,2 2 , Mjf / 2 2 \ A'/f* M/f / 2 , 2 2 \ 



Mb- + mI:- - Mb* mI* + lA- - mI* 



G 

2 ~ r;2 



Mk* 2/i|./i|^.A2 ' 



1 



2Mb* ^3 



G _ mIc* - A*l: 



4 ~ 9 

Mk* 



at; 



rF MB^iMB* / , X , as ai3 + W/iTB* . 2,2 2 2 2 

^1 = A 4(Q;i3+WifB*) + -r -2 2 (Mb* +Mif* -Mb*)(Mk* -Mlf 

(Mk* +Mk)(m|* -Mk*)+mI*(3m^* -Mk)' 



Q!3 2 2^ 



_ "3 MB^_MB|/„2 , „2 „2 ^ A CtlS + ^KB* ^ 



3. Decay BIq B*-f. 



The leading contribution to the coupling constant GssoB'-y comes from the diagram in Fig. 4(a) and is defined by 



G's*oB*7 = 9b'Jb-bk Ibi„bi^ , (A8) 



where Jb* b*-7 is the structure integral 



CXD OO C» 



/ / / rfQ!lrfQ!2rfa3 

^BioBti ^ J J J Af exp(iiij aisif asjcB : 



(A9) 







where 



Bx = -ai2lJ'K - Q^sMs - iA'WbkWkb + " {ij-b^oiibk + MB*a2) • (AlO) 

4. Decay B^i ^ 5^7. 

The leading contribution to the coupling constant Gb^iB,7 is due to the diagram in Fig. 5(a) and is defined by 

GBsiB,^ = ^g^2 jY2 9b,i9b,b*k ^B,iB,7 , (All) 

where Ib^iBsj is the structure integral 

1 V Ai V Md 



oo oo oo / 2 2 

dai(ia2rfQ!3 / 2ai f Md, ~ Mk 





where 



/b.,B37 = / / / J exp(B2) asKB* 4 - 1 + ' ^' ) ) , (A12) 



B2 = -Q!i2Ma: - "sMb* - IJ^b,,wb*kWkb* H r (Mb* aiB*A: + Mb,Q'2) • (A13) 
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5. Decay B*,, ^ B^j. 

The leading contributions to the coupling constant Gs^iB'^y, Fb^iB^-j and HB^iS'-y coming from the diagrams in 
Figs. 6(a) and (b) are defined by 

i?B,iB|7 = Y67r2A2 ^Bai (5k.±^±^ + 5^.0^0 J , (A14) 

where R= G,F ov H and /ij are the structure integrals given by 



oc oo oo 



J J J Af exp(B3) ii , 

exp(B3) L2 , (A15) 





oo oo 'OO 

daida2da3 





OO' oo OO' 



A? 



f f f daida2da3 x r 

J J J A2 e^P(^3) L3 , 

^ 

D 2 2 2 2 I OiSKB' ,2 , 2 \ r A t c\ 

B3 = -OtiilK - a2lJ'K' - "sMb* - IJ'B.i'WB'KWkB' + {iJ-B^iCtiB'K + MB*a2) , (A16) 



and 



L, = 



2+2 A / '^J^/ I I " I 2 A / l^i-l^^ 1 ^ 2 '2a 



Mb* Mb*^i/ Ms* Mb* 

m|. m|* 

L3 = -/J,(i-/3i2), /3i = ^^1^, 02 = ^, /3i2 = /3i+/32. (A17) 

6. Decay BJi ^ S^oT- 

The contribution to the coupling constant Gb^iB'^j, coming from the diagrams in Figs.7(a,b) is defined by 

1 

167r2 A2 



Gb^^bi^^ - ,^2.2 9b.^9bi„ (ffB.±B±^ + 9b'Obo^) ^B^^Bi^i , (A18) 



where Ib^iB'^-j is the structure integral given by 



00 00 00 



lB^iB*„f =111 — ^ exp(i?4)(a3 + WB*if + WBif) , 







2 2 2 2 2 

-B4 = -aif^Q, - a2iJ.B - asfJ-K ~ IJ'b^^'^B'kwkb* - Mb* wbkwkb 

, a^ + WB-K + WBK , 2 ,2 \ 

+ (Ms^i^lJCB* + ^J'B^^(^2KB) , 

A = 2 + ai23. (A19) 
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TABLE I: Decay widths of B^q BsTt and Bsi B^n in keV. The range of values for our results is due to the variation of 
Am from 1 to 2 GeV. 



Approach 


r(B:o ^ Bs^') 


Approach 


r(B,i ^ BtTv") 


Ref. [16] 


21.5 


Ref. [16] 


21.5 


Ref. [20] 


1.54 


Ref. [21J 


10.36 


Ref. [27] 


6.8 - 30.7 


Ref. m 


5.3 - 20.7 


Our results 


55.2 - 89.9 


Our results 


57.0 - 94.0 



TABLE II: Decay widths of B*q —> B*y and Bsi Bsy in keV. The range of values for our results is due to the variation of 
Am from 1 to 2 GeV. 



Approach 


r{B:, ^ b:7) 


Approach 


r(B,i -> Bs-y) 


Ref. 


58.3 


Ref. [16] 


39.1 


Ref. [2^ 


171.4 


Ref. [23] 


106.5 


Ref. [23] 


31.9 


Ref. [23] 


60.7 


Ref. [28] 


1.3 - 13.6 


Ref. [28] 


3.2 - 15.8 


Our results 


3.07 - 4.06 


Our results 


2.01 - 2.67 
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(a) (b) 
FIG. 1; Mass operators of B^g and Bsi mesons. 




FIG. 3: Diagrams contributing to the strong transition Bsi + 7r°. 




FIG. 7: Diagrams contributing to the radiative transition Bsi 
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